Dynamical mean-field theory of correlated hopping: A rigorous local approach 
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A general approach for the description of correlated hopping in infinite dimensions, which is based 
on an expansion over electron hopping around the atomic limit, is developed. Such an approach keeps 
the dynamical mean-field theory local ideology and allows one to calculate the thermodynamical 
functions. A grand-canonical potential functional and a $-derivatible theory that does not introduce 
the self-energy is proposed. As limiting cases the Falicov-Kimball model with correlated hopping 
and the model with broken bonds ("diluted" conductor) are investigated, and the connection with 
the Blackman-Esterling-Berk coherent potential approximation approach in the theory of the binary 
alloy with off-diagonal disorder is considered. 
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I. INTRODUCTION 

A microscopic origin of the different effects and phe- 
nomena specific to transition and rare-earth metals and 
oxides, including high-T c , superconductors is due to the 
strong local correlations between electrons. The sim- 
plest interaction that introduces such local correlations 
is the Hubbard single-site Coulomb repulsion Un^riii. 1 
In addition to such local terms there are others involving 
neighboring lattice sites that describe the charge-charge 
interaction VniUj and correlated or density-dependent 
hopping: 



E 



Hj a h a J° { n is+n jS ) . 



(1.1) 



In general, an effective model with correlated hopping 
can be obtained after integrating out additional degrees 
of freedom 2,3 or as a phenomenological model. It was 
first proposed to describe the properties of mixed va- 
lence systems, but after the discovery of high-T c su- 
perconductivity the interest in correlated hopping as a 
possible new mechanism for a superconducting instabil- 
ity and temperature-induced metal-insulator transition 
increased. 3,4,5,6 In addition correlated hopping is consid- 
ered an important factor in the stabilization of the ferro- 
magnetism and localization of electrons. 7 

Models with correlated hopping, like other models with 
strong electron correlations, cannot be solved exactly. 
However in the last decade, great progress in the descrip- 
tion of strongly correlated electron systems was achieved 
with the development of the dynamical mean-field theory 
(DMFT). 8,9 

The main idea of the dynamical mean-field theory, 
which is exact in the limit of infinite dimensions, is in 
the local (single-site) nature of the self-energy. 9 But this 
statement, as shown by Schiller, 10 is violated for systems 
with correlated hopping, when the self-energy becomes 
nonlocal that breaks the general DMFT scheme. Now 
the effect of correlated hopping is similar to the Hartree 



(mean-field) renormalization of the band energy, as an 
additional nonlocal contribution to the self-energy. 

A possible solution of this problem is in a reformula- 
tion of the DMFT approach without introducing the self- 
energy. It should be noted that self-energy appears in the 
Dyson equation for the one-electron Green's function as 
a correction to the band energy due to many-electron in- 
teractions, which are treated as perturbations. On the 
other hand, one can start, not from one-electron band 
states, but from many-electron local states, and analyze 
the influence of the intersite hopping ty on the creation 
of the band structure. Such a strong-coupling approach 
corresponds to the perturbation theory over electron hop- 
ping. In this case the one-electron Green's function can 
be a solution to the Larkin equation, 11,12 



irn 

Gk(u>) = Sfe(w) + E k (uj)tkGk(ijj), 



(1.2) 



where Ejj(w) is an irreducible part of the Green's func- 
tion that cannot be divided into parts by cutting one 
hopping line, and is connected to the self-energy by the 
following equation: 



■^fc 1 ^) = w + /i - Efe(cj). 



(1.3) 



In the DMFT approach such an irreducible part (ir- 
reducible cumulant) was intoduced by Metzner, 13 who 
proved that it is local in the D — > oo limit. Below it will 
be shown that this statement is more general than the 
one about the local nature of the self-energy. 

In this paper we present a general approach for a 
description of correlated hopping in infinite dimensions 
which is based on an expansion over electron hopping 
around the atomic limit, 13,14 which keeps the DMFT lo- 
cal ideology, and which allows one to calculate the ther- 
modynamical functions. As a limiting case the Falicov- 
Kimball model with correlated hopping and a model with 
broken bonds ("diluted" conductor) are considered. In 
this case, such a strong coupling approach corresponds 



2 



to the generalized locator coherent potential approxima- 
tion developed by Blackman, Esterling and Berk 15 (BEB 
CPA) for the description of the binary alloy with off- 
diagonal disorder. 



II. EXPANSION AROUND THE ATOMIC LIMIT 

In general, the hopping term of the Hamiltonian with 
correlated hopping for the Falicov-Kimball and Hubbard 
models is written 

H t = -j= E [ tld l d i + t2d\dj (n if + n jf ) 



It is obvious that 



to) 



(2.1) 



+ hd\djn if n j} 



and 



H t = 77= E [*i a L a JT + t24 a a j<7 {n lS + n jd ) 



(2.2) 



respectively. As a rule, the first two terms with t\ 
and t 2 were considered, 4,5 ,7, 10 but an influence of the 
last term with i 3 , that can originate from the nondi- 
rect interactions over other degrees of freedom, was 
also investigated. 3 ' 6 ' 16 In Hamiltonians (2.1) and (2.2) 
the terms t\, t 2l and t 3 describe hopping between two 
nearest-neighbor sites i and j irrespective of the occupa- 
tion of these sites, the hopping between these sites when 
one of them is occupied by another particle and the hop- 
ping between two occupied sites, respectively. One can 
transfer Hamiltonian (2.1) to the other one that describes 
the hopping between the sites with the certain occupa- 
tions, 

H * = "if E [% + P?4djP? + tifPfdldjPr 



(ij) 

t+r P +d\d 3 Pr +tr+pr d \d j pf 



(2.3) 



where the projection operators 

Pf = n if , Pr = l-n if 
and matrix notations 



d, = 



t ++ tfr 
d t — I V , 13 



(2.4) 



(2.5) 



are introduced. For the Hubbard model it corresponds 
to the transition to the Hubbard operators 



Ht = -U a ^ = y^L ) ■ (2-6) 



(2.7) 



and hopping integrals in Eqs. (2.1), (2.2), and (2.5) are 
connected by the following equations 

t— = t u h = t—, 

t+~ =t~+ =h +t 2 , t 2 =t+-(-+^ -i— , (2.8) 
t ++ = ii + 2t 2 + i 3 , h = t ++ + r~ - t+- - t~ + . 

The total Hamiltonian of the electronic system with 
local correlations and correlated hopping, 



H = Y,H l + H t 



(2.9) 



includes, besides the hopping term H t that is not local, 
the single-site contributions Hi. Our aim is to consider 
the D — > oo limit, and, as mentioned in Sec. I, it is 
convenient to start not from the Dyson equation, that 
considers the terms with correlated hopping t 2 and i 3 
as some kind of many-particle interactions, but from the 
Larkin equation 11,12 in a coordinate representation, 

Gij(u) = Hjj(w) +^S i j(w)tj m G m j(w), (2.10) 

Ira 

that treats all hopping terms in the same manner. Here 
an irreducible part (irreducible cumulant 13 ) Hjj(a>), that 
cannot be divided into parts by cutting one hopping line 
t; TO , is introduced. 

For the models with correlated hopping all quantities 
in Eq. (2.10) are matrices, e.g., the components of the 
Green's function Gy(w) are constructed by the projected 
(Hubbard) operators 



Gt^ --') = - (T(P*dMd]P]) T , / - " a-rap _ T y 



sn 



G«(t -t') = - (Td 4 (r) ® dj(r')) = (3 



sn 



(2.11) 



where O is a grand canonical potential functional and 
St J" (t' — r) is some auxiliary field. 

After the Fourier transformation to the momentum 
space one can find the Green's function as a solution of 
Larkin equation (2.10), 



G fc (w) = [1-S fc (w)t fc ] 'HfcH 

= [H fc 1 H-t fc ]- 1 , 

and the total physical Green's function 



G^t-t') 



Td^dUr') 



Shji(T' - t) 



is equal to the sum of all matrix components 



(2.12) 



(2.13) 



(2.14) 
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In general, an irreducible part Sjj(u>) is represented 
diagrammatically by the single-site vertices (because all 
interactions in Hi are local) connected by hopping lines 14 
and in the D — > oo limit it can be shown (see, e.g., 
Ref. 13) that it becomes local 



j(u>) = 5ijS(u), a k (uj) = S(w). 



(2.15) 



Such a matrix representation allows one to reformulate 
the dynamical mean-field theory of systems with corre- 
lated hopping in terms of local quantities. Indeed, the 
local irreducible part S(w) depends on the electron hop- 
ping only via the local coherent potential (auxiliary A 
field of Brandt and Mielsch 17 ) 



or, in matrix notations, 



(2.16) 



(2.17) 



where Gj^(w) is the Green's function for lattice with 
the removed site o. On the other hand, we can introduce 
another quantity 

I(W) = G '™M t mo = -J- G fc(^) *fe ( 2 - 18 ) 



Im 



that is connected with coherent potential [Eq. (2.17)] by 
the equation 



I(w) = [1 - JHSfwjp 1 J(w). 



(2.19) 



From Eqs. (2.17), (2.18) and (2.19) we obtain the follow- 
ing equation for the coherent potential J(u>): 

1 ^ [H- 1 ^) - t fc ] - 1 = [H- 1 ^) - J( W )] 1 , (2.20) 

that is a matrix generalization of the known one (see, 
e.g., Ref. 9). 

The expression in the right-hand side of Eq. (2.20) is 
the Larkin representation of the Green's function, 



G imp H=[S- 1 (w)-J(w) 



(2.21) 



for the effective single-impurity problem with the statis- 
tical operator 



Pimp 



- pHe Texpl-JdT JdT'dt{T)J(T -r')d (T')> • 

loo J 

(2.22) 

Equations (2.20) and (2.21) give the closed set of equa- 
tions for the local coherent potential J(w) and irreducible 
part H(w), 

1 ]T [G im » + J(w) - t fc ] 1 = G imp (c), (2.23) 



(2.24) 



that was derived without introducing of the self-energy 
for the total Green's function [Eq. (2.14)]. Solution 
(2.24) is substituted into the expression for the lattice 
Green's function [Eq. (2.12)] that allows one to calculate 
the total Green's function (2.14). Below it will be shown 
for the Falicov-Kimball model with correlated hopping 
how to solve the effective single-impurity problem, to cal- 
culate Gi m p(w), and to derive the self-energy that is now 
momentum dependent. 

Another advantage of such an approach is in the cor- 
rect analytic properties of all quantities. In the complex 
plane they behave as 



(2.25) 



for \z | — > oo, that allows one to build the Lehmann rep- 
resentation for the Green's functions 



+ oo 



i \ 1 f j ImG imp (w - i0+) 
G imp (z) = - / doj- 



z — UJ 



+ oo 

1 f , ImG fc (^-iQ+) 
Gfc(z) = - / duo . 



z — uj 



(2.26) 
(2.27) 



the irreducible part 



3(z) 



+ 00 

— duj 

71" J 



z — UJ 



(2.28) 



and the coherent potential 



+ oo 

J(*)=i /^ ImJ( - i0+) . 

7T ,/ Z — UJ 



(2.29) 



The last two equations follows from the definitions of the 
irreducible part [Eq. (2.10)] and the coherent potential 
[Eq. (2.17)]. It should be noted that we can, if we wish, 
introduce a matrix of the local (CPA) self-energies 15 in 
the same way as in Eq. (1.3) by 

S(z) = (z + /i)l - H _1 (z), (2.30) 

but components of such self-energy matrix do not possess 
correct analytic properties because C£ a ^ 1 and they 
diverge for \z\ — ► oo (see below for the Falicov-Kimball 
model). 



III. SUMMATION OVER WAVE VECTOR 

In Eq. (2.20) the sum over wave vector k can be cal- 
culated directly when we put 



(3.1) 



4 



where a is a matrix of the correlated hopping constants. 
In this case, in (2.20), we can replace the sum over k by 
the integral with some density of states p(t) 



+00 



h 

K -00 
+ 00 

= J dtp(t) [a _1 E _1 (w) -il] _1 a _1 . (3.2) 

—00 

Matrix M = a _1 H _1 (o)) can be diagonalized by some 
transformation V, 



V^MV : 
and finally we obtain 



Zl 








Zs 



(3.3) 



G imp a = V 







J = aV 









F(zs) 




Zs - 



v- 1 , 



v 



(3.4) 



(3.5) 



where 



+00 



F(z) = / dt 



z-t 



(3.6) 



is the Hilbert transform of the density of states p(t). 

Now, let us consider some popular densities of states. 
The Hilbert transform of the Lorentzian density of states 



is equal to 



p(t) 



F(z 



W 1_ 

~K W 2 + t 2 



1 



' ztiW 
that immediately gives 

J(lu) = ±iWa. 



(3.7) 



(3.8) 



(3.9) 



For a Bethe lattice with a semielliptic density of 
states, 9 



(3.10) 



we have 



^) = ^2 (z-V^w*), (3.ii) 

(3.12) 



1 w 2 

m - - F{z) 



that gives a matrix generalization of the known relation 18 

(3.13) 



W 2 

J(w) = — aG imp (w)a. 



It should be noted that Eqs. (3.4) and (3.5) are correct 
only when deta 7^ 0. In other case, a direct decompo- 
sition of Eq. (3.2) into simple fractions over tk must be 
used. But results for the particular cases [Eqs. (3.9) and 
(3.13)] are always correct. 

For a density of states different than the Lorentzian 
and semielliptic ones the iterative procedure is as fol- 
lows. We start from the initial value for the coherent 
potential 3(u>), e.g., that of the Lorentzian density of 
states [Eq. (3.9)], that is used to calculate the Green's 
function Gi mp (o>) for the single impurity problem [Eq. 

(2.22) ] (see below for the Falicov-Kimball model). Then 
Eq. (2.24) is used to find an irreducible part E(w), Eq. 

(2.23) is employed to find a new value of Gi mp (w), and 
then Eq. (2.24) is used to find a new coherent poten- 
tial J(oj) value. The algorithm is then repeated until its 
converges. 



IV. THERMODYNAMICS AND 
cfr-DERIVATIBLE THEORY 

One can see that the DMFT system of equations pre- 
sented above is a matrix generalization of the known 
equations for the systems without correlated hopping, 
with the replacement of the local self-energy by the local 
irreducible part 



S(w) = uj + p-^{uj) 



(4.1) 



However, there are principal difficulties in obtaining ex- 
pression for the grand canonical potential (free energy). 

The main idea in derivation of a grand canonical po- 
tential in the DMFT is the following. 9 ' 17 We start from 
the Baym-Kadanoff functionals for the lattice, 



N 



v k 

~]]Y.-^J2 G k {iu, v )E(iu) v ) + 

v k 

and the effective single-impurity problem: 

^im P = - ^2 In [iuj v + p - ^(ioJ v ) - J{ioJ v )\ 



(4.2) 



4 ^2 Gi mv {iu v )Y,{iu u ) + $ imp . 



(4.3) 







In the case of the absence of correlated hopping the self- 
energies for the lattice and impurity are local and the 
same. In addition 



(4.4) 



■5 



and 



-y 



G k {iui v ) = G imp (iu> v ), 



(4.5) 



that allows one to exclude an unknown Luttinger-Ward 
functional and to find the final expression 



= fiimp - 4 e{ ji E ln \ iujv + v ~ - 

i> ^ k 

+ In G imp (iuj v )\. (4.6) 



But in the case of correlated hopping the self-energy is 
momentum dependent 10 and it is impossible to write 
down the Baym-Kadanoff functional for the impurity in 
the form of Eq. (4.3), and that is the main problem. 

In order to solve this problem let us construct the 
Baym-Kadanoff-typc functional without introducing the 
self-energy To do this we start from the expression for 
the Green's function (2.12) and (2.11) in the Larkin rep- 
resentation 



l|GT(« 



lat 



[a-s fe (w)tfe] 1 s k (u), 



(4.7) 

where we introduce an auxiliary field t k (u>) = t k +6t k (w), 
and a known relation 



S ln det A = Tr (A' 1 6 A) . 



(4.8) 



That gives 



^lat „ , ^' 



(4.9) 



7? E 7j E lndct t 1 ~ &k{iuv)tk(iuv)] ■ 



v k 



where Oo is the grand canonical potential for the single- 
site Hamiltonian Hi (tjj = 0) and 



a f gfc^oy) 



x tfe- (iuv ) [1 - Sfe- (iuv )t fe < (iuv )] 1 1 . 



(4.10) 



As mentioned above, the irreducible part Hfc(iu„) can be 
diagrammatically represented by the single-site vertices 
connected by hopping lines t k {iuo u ). 14: On the other hand, 
by summing up the series of the diagrams of the same 
topology, one can see that hopping lines are collected 
into chains, 

tfe («*;„) = t k {iuj v ) + t k {iuj v )'S k {iuj v )t k {iuj v ) + ■■■ 
= t k (iu u ) [1 - B, k (iuj v )t k {iijO v )}~ 1 , 

(4.11) 



and the irreducible part can be represented by some 
skeletal diagrams with the t k (iu>„) replaced by t k (iu> v ) 
(spreading of lines) . Now we can represent O' as 

| = -^E^E Tr [b*(^)«*M + (4-12) 



where $i at is a Luttinger- Ward-type functional: 

<5$lat 



07? 



(4.13) 



Finally, for the grand canonical potential functional for 
the lattice we obtain 



= fi 0~ -^^2\ndet[t - a k (iuj u )t k (iuj v )} 

P v k 



(4.14) 



In the limit of infinite dimensions an irreducible part 
becomes local [Eq. (2.15)], and from Eq. (2.20) it follows 
that 



-^^t fc (^) = 3{iw v ), 



where 



3(iu v ) = 3(iw v ) [1 - a{iuj v )3{iuj v )] 1 . 



(4.15) 



(4.16) 



On the other hand, we can write a functional for the 
impurity in the same form, 



ft, 



=f?o — 4 E m< ^ ct — H(zcj„)J(«cj„) 

-^E Tr [ H( ^ )5( ^ } 



(4.17) 



where 



<5JT Q (i^) 



S^(z^). 



(4.18) 



From Eqs. (2.15), (4.13), (4.18), and (4.15) it is easy to 
prove that 



flat 

N 



imp 7 



(4.19) 



and finally we obtain an expression for the grand canon- 
ical potential for the lattice in terms of the quantities for 
the impurity problem 



filat 
~N~ 



1 f 1 

= ^imp ~ Tj Ys\ N E llld<3t t 1 ~ H (^) t fe] 

— lndct [1 — E(iu v )3(iu) v )] 



(4.20) 
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In the absence of correlated hopping this reduces to the 
known expression [Eq. (4.6)]. 

Logarithms in Eq. (4.20) originate from the sum of the 
loop diagrams, 



and 



Tr^ y [3(iuj v )J(iu} u )] 1 



(4.21) 



(4.22) 



and using analytic properties of the irreducible part [Eq. 
(2.28)] and the coherent potential [Eq. (2.29)], in Eq. 
(4.20) one can replace the sum over Matsubara's frequen- 
cies by an integral over real axis, that gives 



''imp 



7T / + 1 



x ji^Imlndet [l - S(w - »0 + )t fc ] 

^ k 

- Imlndet [l - S(w - iO+)J(uj - i0+)] J. 

Expressions (4.14) and (4.13) are very similar to the 
one known in the Baym-Kadanoff approach 19 but now, 
instead of the conjugated quantities the Green's function 
Gk{iu v ) and self-energy Jlk^iuj^), we have the renormal- 
ized hopping tk(iu) v ) and the irreducible part 5k(iuj v ). 
Such a formal analogy allows us to build a <I>-derivatible 
theory starting from functional (4.14) in the same way 
as done by Baym and Kadanoff. 

First of all, for the electron concentration we have 



Hi 



1 dfl 



Iat 



imp 



N d/ii d/j,i 

cMp 1 gglat 

dfii N dfii 1 



J— const 



(4.24) 



where a partial derivative is taken over /i not in chains 
tfcCiw,) [Eq. (4.11)]: 



d _ d 
dfi d/i 



ajuk 



dt1 1 {iuj y ) 5 
dp 5^{iu) v )' 



(4.25) 



Next, for the correlation function (dynamical suscep- 
tibility) 



X12 



dri\ dri2 
d(j> 2 d[ii 



i d 2 n 



lat 



N d\i\d[ii 



(4.26) 



we obtain an expression that can be represented diagram- 
matically as 20 



Xl2 



a -'0 2+, CD3 2 ' <427) 



where the full four-vertex is a solution of the Bethe- 
Salpeter type equation 



□ -□-LTD 



(4.28) 



Here the thick wavy lines represent the renormalized hop- 
ping tk(iu) v ) [Eq. (4.11)], and 



d 2 n i a 2 $ lat dn x dn 2 



d (5$i at 8Z 23 Sni 



1 <1 =P— 

N 3 dpi St 32 9fii St 32 

□ 



(4.29) 



1 I 14 

8t 2 i8t i3 



<5 2 $i at 5E U 5S, 



■34 



5*21 



are irreducible vertices that can not be divided into parts 
by cutting two hopping lines, and it can be shown that in 
the D — ► oo limit all irreducible vertices become local. 20 



V. FALICOV— KIMBALL MODEL 

Now let us apply the developed above approach to the 
Falicov-Kimball model with correlated hopping. In this 
case the single-site Hamiltonian Hi is written as 

Hi = -fifUif - fidriid + Un ld n if (5.1) 



or 



Hi = [(U - Hd)nid - M/] p ? - fJ-dfiidPi 

= H+P+ + Hr P r. ( 5 .2) 



Projection operators and P~ commute with the total 
Hamiltonian [Eq. (2.9)] and the statistical operator for 
the effective single-impurity model [Eq. (2.22)] can be 
represented as 



Pimp = P+p++P p_, 



(5.3) 



where 



p a = e-^Texp i-J dr J dr> 4(r) J aa (r - r')d (r')| 

(5.4) 

are statistical operators for the noninteracting d electrons 
placed in the different fields for the different subspaces 
a = ±. As a result, the partition function for the impu- 
rity is a sum of the partition functions for the subspaces: 



■^imp Tr Pimp & 



1 



_ e -(3Q+ , e -PQ- 



(5.5) 



Q+ = - t if - - In [1 + e-l 3 ^-^ 
l Eln / J++(^) 







icv v + [id - U 



(5.6) 



Q_=-iln(l + e^)-i£ln(l 
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that gives, for the grand canonical potential of the impu- 
rity, the following expression 



Here 



Iln(e-W++e- w -). 







(5.7) 



n ^ = ^ P ) = e -/3Q+ + e -/3Q- 



(5.9) 



n ^ ._ „. . , . , , , is the /-particle concentration. In addition, for the d- 

hom fiq. (5.3) it also follows that components ol the ,. , ... , 

. r \ ' r ,, . ., K particle concentration we have 



Green's function for the impurity are equal 

(P+) 



G impM = 



u + Hd-U - J++{uj) 

(P-) 
10 + [Id - J (to) ' 



++ 

imp 



{iu v ) + G-Uu v ) . (5.10) 



(5.8) 



Gt~M = Gr+ (w) = 0. 



The coherent potential J(w) is a solution of Eq. (2.23), 
and finally for the Green's function of the lattice we ob- 
tain 



G fc (w) 



V k (uj) 



[u> + H d - J-(w)(P+) t k -{p-)] (P+) [J+-( W ) - *+"] (P+)<P-) 

- [J-+M - t fc +] (p+)(p-) [ u + p d -u- j++(w)(p-) t++(p+)] (P-) 



(5.11) 



v h (u) = [u + nd- j (uj)(p + ) - t k -(p-)] [ W + Iid -U- j ++ (cj)(p-) - 4+(P + )] 

The inverse matrix of the irreducible part is equal, 



^[w + ^-tf- J++( U )(P-)] 



J"+( W ) 



J+-( W ) 
i + J— (w)<P+)] 



(p-> 



(5.12) 



(5.13) 



and one can see that in the \u\ — > oo limit its diagonal 
components diverge faster than w that does not allow to 
introduce a matrix of the local (CPA) self-energies with 
correct analytic properties by Eq. (2.30). As a rule, 
the behavior of the CPA self-energy at the u> — > ±oo 
limit is not analyzed considering analytic properties of 
the BEB CPA. 21,22 On the other hand, the components 
of the irreducible part matrix at — > oo behave as 



(P H 



3+-(c,).E- M 



(P- 



U) 



(5.14) 



that is in agreement with Eq. (2.25). 

In the absence of / particles (P + ) = and (P~) = 1, 
and we obtain solution for the free electrons, 



(5.15) 



whereas for the opposite case, when all sites are occupied 
by /-particles (P+) = 1 and (P") = 0, 



1 



u> + ii d - U - 4 + ' 



(5.16) 



In general the total Green's function [Eq. (2.14)] for 
a Falicov-Kimball model with correlated hopping can be 
written in the Dyson representation as 

GfcM = — 1 n r- . (5-17) 



w + - E fe (w) - ife' 



where 



(5.18) 

is the Hartree (mean-field) renormalized hopping already 
introduced by Schiller, 10 and 

£ fc (w) =U(P+) + J 3 (uj)(P + )(P-) (5.19) 
Uik(u)U 2k (u;){P + )(P-) 

LO + fM - U(P-) - J(UJ) - t 3k (P+)(P-) 

is a momentum-dependent (nonlocal) self-energy. Here 

CM") =^ + K + - t k + - J—(«>) + J~ + H] (p + ) 
+ [tt--t k --J + -n + j++(c)] (P-) 

(5.20) 

U 2k (u) =U + [t++ - t+- - J (w) + J+~M] (P+) 
+ [^ + -^--J-+(o;) + J ++ M] (P-) 
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describe renormalization of the interaction U, and 
J(w) =J ++ (cu)(P-) 2 + J—(uj)(P+) 2 (5.21) 

+ + (p+xp-), 

J 3 (lo) =J ++ (lu) + J~~{uj) - - J~ + (u), (5.22) 



*3fc -tt + +t k - 4" 



t fc + (5.23) 



In his paper Schiller 10 considered the case of 

hk = 0, (5.24) 

and from [Eq. (5.19)] we obtain the same result for the 
self-energy, 



E fc (w) = S (w) + Si(w)t fc + S 2 (w)*fc. 



(5.25) 



that for the nearest-neighbor hopping i*. contains only 
local, nearest-neighbor and next-nearest-neighbor contri- 
butions. The connection with Schiller's approach for this 
case is given in the Appendix. But in a general case of 
hk 7^ 0, when Schiller's approach is not applied, the self- 
energy is extended over the whole lattice, see Ref. 22 for 
a comparison with the BEB CPA. 

In Fig. 1 we present the spectral function (imaginary 




FIG. 1: The spectral function (imaginary part of the Green's 
function) and its components for different relations ti /ti {tz = 
0) for the D — co hypercubic lattice with nearest-neighbor 
hopping (W = 1, U = 2, T = 0.01) at half-filling fi d = fj, f , 
rif + n d = 1. 

part of the Green's function) and its components for dif- 
ferent relations t 2 / 1\ (t 3 = 0) for the Gaussian density of 



states (D = oo hypercubic lattice with nearest-neighbor 
hopping), 



P(t) = 



1 



W^/tt 



(5.26) 



(W = 1) at half- filling [i^ = /x/, n/ + rid = 1- One can 
see that these results are in good agreement with the one 
obtained by Schiller 10 for the D = 2 square lattice. At 
points t2/t\ = and £2/^1 = -1 we have a symmetric 
band picture, whereas for £2/^1 > and ^Ai < —1 pa- 
rameter values of the correlated hopping the upper band 
become wider and for — 1 < t 2 /ti < values the lower 
band does so. Such a behavior can be explained in the 
following way. The symmetry of the band picture is de- 
termined by the relation \t ++ \/\t | and by the deviation 
of the /-state occupation from half-filling rif — ^ . For the 
symmetric points we have t ++ = t~~ for the £2/^1 = 
case and t ++ = —t~~ for the ^Ai = — 1 case. In both 
cases rif = i for any temperature value, but out of these 
two points there is a deviation from the half-filling of 
the /-state occupation (see Fig. 2), and this deviation 



0,10 




tA 



FIG. 2: Deviation from the half-filling of the /-state occupa- 
tion as a function of t2/ti. Parameter values are the same as 
in Fig. 1. 

achieves its maximum at low temperatures. On the other 
hand, at high temperature n/ — ► | for T — ► 00. The gap 
in the spectrum is at a maximum at rif w ^, and with 
the decrease of the temperature such a deviation from 
the half-filling of the / states can lead to a shutdown of 
the gap and a Mott transition (see Fig. 3). 

Now let us consider some limiting cases. In the absence 
of correlated hopping (i 2 = £3 = 0), 



t 



++ 



tl 



t k = t k 



(5.27) 



and 



J++H = j—(u) = J + - 



(u) = J-+(w) = J(w), 

(5.28) 
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T=0.01 

T=0.1 

T=1 




FIG. 3: Temperature development of the gap (Mott transi- 
tion) for tz/ti = 0.5, U = 1.7 (top) and t 2 /h = -0.5, U = 0.3 
(bottom). 



that follows from definition (2.17), and we obtain the 
known result for the Falicov-Kimball model, 



E(w) = U(P+) + 



tk — tk, 

U 2 (P+)(P- 



< + H d -U{P-)-J(u) 



(5.29) 
(5.30) 



with local self-energy. 

Another symmetric point t 2 /ti = — 1 corresponds to 
the case of diagonal hopping matrix (t^~ = t k + = 0) 
when hopping is allowed only between sites with the 
same occupancy. In this case we have a binary alloy with 
atoms of two types distributed stochastically over the lat- 
tice, and hopping is restricted between atoms of different 
types. Now atoms of one subsystem can be treated as 
impurities for the other subsystem, and it was supposed 
in Ref. 15 that this case corresponds to two independent 
bands (independent band limit). In this limit the particu- 
lar case of ts = was considered for the one-dimensional 
Hubbard model when exact results can be obtained. 23 
In Ref. 24 it was shown for the Hubbard model, by the 
slave-boson mean-field approach, that in this case a direct 



transition from the superconducting state to the Mott in- 
sulator takes place. 

For the Falicov-Kimball model and the semielliptic 
density of states [Eq. (3.10)], a coherent potential matrix 
in this case is also diagonal [J -1 (w) = J h (w) = 0], 



J—(w) = 



(5.31) 



4 

x - uj + vd-u± %yjwi - (a, + Md - uy , 

a 2 W 2 

— — G—(w) (5.32) 



u + f i d ±tJW 2 - (uj + ^ d ) 2 



w ± = a ±± w2 ( p± )- 

The Green's function for the lattice contains two contri- 
butions 



G + k + = 
G7- = 



(P+) 



cu + fM -U- J++(u)(P-) - t k (P+) ■ 

(PJ 

w + \ia - J—(u){P+) - tk(P-) 



(5.33) 
(5.34) 



that describe two independent bands originating from 
two disordered subsystems characterized by their own co- 
herent potentials. The spectral weight function contains 
two bands 

p(w) =- Im G imp (oj - ii - i0 + ) 
2 



separated by the gap 

A = U- \W+\ 



\W-\, 



(5.35) 



(5.36) 



that, in general, is temperature dependent and disap- 
pears at some critical value 

U c = \W+\ + \W-\ = W (|a ++ |v^) + |a__|v^P=)) 

(5.37) 

(the Mott transition point) that achieves its maximum 
CI,™ = Wy/a? ++ +a?__ (5.38) 

at 



a\ + + a 2 _ 



(5.39) 



At half-filling (/if — fid, nf + n d = 1) and for the case 
^3 = 0, when a__ = 1, a ++ = —1, and nf = rid = i, the 
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gap does not depend on temperature and disappears at 

u c = wVz- 

It should be noted that this result for independent 
bands is correct only for the scmiclliptic density of states 
when coherent potential 3(lj) and diagonal Green's func- 
tion for the impurity Gi mp (w) are connected by relation 
(3.13). But in a general case it is not correct, e.g., for 
the Gaussian density of states when each subband con- 
tains contributions from both subspaces (see Fig. 1, case 

t 2 /h = -l). 

The next simplification of the model corresponds to a 
system with broken bonds when we have hopping only 
between the sites not occupied by / particles 



t ++ - t + ~ - r 

l k ~ L k — k 



= 0, 



J++H = J + -(w) = J _+ (w) = o. 



(5.40) 
(5.41) 



In other words, lattice sites are occupied stochastically by 
atoms of two types ( "conducting" and "isolating" ) , and 
electron hopping is allowed only between atoms of the 
first type and the influence of other subsystem is consid- 
ered as disorder (the case of the "diluted" conductor). 
Now the total Green's function [Eq. (5.17)] contains the 
renormalizcd hopping 



tk — t. 



(P- 



and some momentum-dependent (nonlocal) self-energy 
E fc (w) = U(P+) + J—(uj)(P + )(P-) (5.42) 

[u j- - h (p+) -t k -(p-)} 2 (p+) (P-) 

■ + W - U(P-) J-{u>)(P+)* - t^-(P+)(P-) ■ 



+ 



But, in fact, these quantities have no any physical sense. 
Indeed, in this case the total Green's function contains 
two contributions 



u + ii d -U 



(P- 



uj + fx d - J—(w)(P+) - t k - (P-) 



(5.43) 
(5.44) 

(5.45) 



where G k + {to) originates from the localized levels of the 
isolating atoms and G k ~{uo) describes the one-electron 
excitations in the disordered conducting system charac- 
terized by the coherent potential J~~(uj). For the semiel- 
liptic density of states [Eq. (3.10)], we have 



and the one-particle spectral weight function 
1 



(5.46) 



■ ImG imp (w — ft — 



(5.47) 



= (p+)6(w -u) + (r)-^/i^ 



contains the 5-peek from the localized levels and con- 
ducting band of the half- width W- — W \J (P~) and 
their statistical weights are equal to the concentrations 
of the isolating (P + ) =nj and conducting (P~) = 1 — n/ 
atoms, respectively. A critical value of U when a "gap" 
develops is equal, 



(5.48) 



and is a function of the temperature. This limiting case 
of the Falicov-Kimball model with correlated hopping 
shows that the explanation of the model behavior in 
terms of the CPA approach is more appropriate than in 
the Fermi-liquid language. 



VI. SUMMARY 

In this paper we presented a general approach to the 
description of correlated hopping in dynamical mean-field 
theory, whose application is not limited by the case of 
the Falicov-Kimball model and can be applied to other 
models, e.g., Hubbard, t-J, etc. It is based on the Larkin 
equation (an expansion over electron hopping around the 
atomic limit) that considers all hopping terms, includ- 
ing correlated hopping, in the same manner. Another 
starting point is the local character of the irreducible 
part (irreducible cumulant) of the Green's functions, con- 
structed by the projected (Hubbard) operators, in the 
D — > oo limit, that is a more general statement than the 
one about the local character of the self-energy, which 
in the case of correlated hopping is nonlocal. Such an 
approach keeps the dynamical mean-field theory local 
ideology, and allows one to calculate the thermodynam- 
ical functions. To do this a grand canonical potential 
functional of the Baym-Kadanoff-type, that allows one 
to build a $-derivatible theory without introducing the 
self-energy, was proposed. 

As an example, a Falicov-Kimball model with corre- 
lated hopping, and its limiting case of a model with bro- 
ken bonds ("diluted" conductor), were considered when 
exact analytic results can be obtained. In particular, 
the temperature-driven Mott transition in the Falicov- 
Kimball model with correlated hopping was investigated. 
Such a strong coupling approach for the Falicov-Kimball 
model with correlated hopping corresponds to the BEB 
CPA approach for a binary alloy with off-diagonal disor- 
der, and it follows that the BEB CPA becomes exact in 
the D = co limit. In addition, analytic properties of the 
local self-energies in the BEB theory were considered. 
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APPENDIX: EQUIVALENCE TO SCHILLER'S 
APPROACH 



A solution of this system of equations with the use of Eq. 
(A. 2) gives 



For the case t^k = we can find a connection between 
components of £ and J (here and below we use J H = 
J h and T,^ d = by comparison of Eqs. (5.25) and 
(S.8) [here and below (S.8) indicates Eq. (8) in the paper 
by Schiller 10 ]: 



J++ 



U+ E 



Md 



dd 



u 



u 



i^d-^r 



yS 



yS 



J+- = 



yS 



yS 



n f 



l-n f ' 
(A.l) 



In addition, according to their definition the components 
of the Green's function (S.9) are not independent and 
satisfy the following relation: 

( w - ^dd)Gdd - (h + 2S^ d )G^d - Y.^Gw = 1. (A. 2) 

On the other hand, from Eqs. (2.20), (2.21), (5.11), and 
(A.l), we obtain a connection between the components 
of the Green's functions 



G++ = 



yS 



/ 



n f 



^4,d 



yS 

^^ptp 



Gdd — hG^d 



y s j_ (H d -t 2 y ' 



U 



G = 



yS 



+ 



Vld-h? 



(ti + 2ti)G^d 



yS 



G d d 



(^) 2 



~ ^dd 



(A.3) 



G+- = 



yS 



(ti + t 2 )G^ d 



J dd 



yS 

f^t±/yS 

-y-s-v^ipd 



(^W - *2) ] Gdd 



0. 



G dd = G++ + G—, 
G^,d 1 



^4>d ~ t2 G ++ _ Z2t±Q~ 



yS 



yS 



(A.4) 



G 



tplft 



yS _ tj yS 
^dd u Q++ , n dd Q- 
yS 

^^pip ^iitp 



that is equivalent to the solution of Eqs. (S.5) and (S.10) 
with respect to G: 



G = -n f (± s ) '-(l-nf) 



U t 2 
t 2 



. (A.5) 



For the /-particle occupation number rif = l/(e^ £ + l), 
where e is defined by Eq. (S.ll), we have 



e-E f + fi- 



ll 



rij det 



U t 2 

t 2 o 



t S {lLO v ) 



(l-n/) 2 detE s (^) 
iuo v + fi — U — J ++ (iuo„) 
iuj v + fi — J (iw„) 



(A.6) 



On the other hand, from Eqs. (5.9) and (5.6) we have 



1 1 + e~ /3 ( c/ ~' id ) 
-Q_ = -^--ln (A.7) 



^5> 



(3 1 + e^" 
IL0 V + /J, d — U — J++ (iuj v ) lLO v + fid 
iuj v + fid - J iu v + [id - U 



that gives the same result if we use 



lVln^ + ^-^lln 1 + e '7 j) . (A.8) 
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